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Abstract 

In this paper, novel closed-form expressions for the level crossing rate (LCR) and 
average fade duration (AFD) of k — fx shadowed fading channels are derived. The new 
equations provide the capability of modeling the correlation between the time derivative 
of the shadowed dominant and multipath components of the k — /j shadowed fading 
envelope. Verification of the new equations is performed by reduction to a number of 
known special cases. It is shown that as the shadowing of the resultant dominant 
component decreases, the signal crosses lower threshold levels at a reduced rate. 
Furthermore, the impact of increasing correlation between the slope of the shadowed 
dominant and multipath components similarly acts to reduce crossings at lower signal 
levels. The new expressions for the second-order statistics are also compared with field 
measurements obtained for cellular device-to-device and body centric communications 
channels which are known to be susceptible to shadowed fading. 
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I. Introduction 

The k — [i shadowed fading model first appeared in the literature in [Ql and immediately after 
this in fl2|. It has been proposed as a generalization of the popular n — [i fading model 0. In this 
model clusters of multipath waves are assumed to have scattered waves with identical powers, 
alongside the presence of elective dominant signal components - a scenario which is identical to 
that observed in k — n fading OJ. The key difference between the k — fi shadowed fading model 
and that of classical n —/i fading is that the dominant components of all the clusters can randomly 
fluctuate because of shadowing. In particular it is assumed that the shadowing fluctuation follows 
a Nakagami distribution 01. Like the k — // distribution, the k — // shadowed distribution 
is an extremely versatile fading model which also contains as special cases other important 
distributions such as the One-Sided Gaussian, Rice (Nakagami-n), Nakagami-m and Rayleigh 
distributions. In addition to this, it also contains as a special case Abdi’s signal reception model 
||5l which considers Ricean fading where the dominant component also undergoes shadowed 
fading which follows the Nakagami distribution. Due to the ability of the Nakagami probability 
density function (PDF) to approximate the lognormal PDF [j6]|, the k — // shadowed fading model 
can also be used to estimate Loo’s well-known model for land mobile satellite communications 

d 

While the research of composite fading models such as the k — [i t gamma model [[8J and 
its associated second-order statistics including the level crossing rate (LCR) and average fade 
duration (AFD) J9l have been advanced, unfortunately, at present, similar closed-form expressions 
for the second-order statistics of the k — /r shadowed model are currently unavailable in the 
literature. The LCR and AFD of a fading signal are of great importance in the design of 
mobile radio systems and in the analysis of their performance iflOl . Among their many potential 
applications are the design of error correcting codes, optimization of interleaver size and system 
throughput analysis as well as channel modeling and simulation. In this paper, convenient closed- 
form expressions for the level crossing rate and average fade duration of k — [i shadowed fading 
channels are derived and subsequently verified by reduction to known special cases. An important 
empirical validation is also performed through comparison with field measurements from two 
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different types of wireless channel which are known to suffer from shadowed fading, namely 
cellular device-to-device (D2D) communications channels [j2] and body centric communications 
channels m. 

The remainder of this paper is organized as follows. Section II provides a brief overview of 
the k — /i shadowed fading model. Important relationships between the k — /i shadowed fading 
envelope and its time derivative, which underlie the second-order equations proposed here, are 
established in Section III. Also presented in Section III is the derivation of the LCR, while the 
derivation of the AFD is given in Section IV. The new expressions for the LCR and AFD of 
the k — /i shadowed fading model are compared with some empirical data obtained from field 
measurements in Section V. Lastly, Section VI finishes the paper with some concluding remarks. 

II. An Overview of the k - /i Shadowed Fading Model 

The k — fi shadowed fading model was originally proposed in [l]. A slight variant of the 
underlying signal model was also developed independently and appeared shortly after this in 
i. In m a rigorous mathematical development of the model was performed, while in [O the 
model was the result of channel measurements conducted to characterize the shadowed fading 
observed in device-to-device communications channels. Both papers have developed important 
statistics related to the k — ji shadowed fading model such as the probability density function 
and moment generating function. In [[0, the cumulative distribution function (CDF) and the sum 
and maximum distributions of independent but arbitrarily distributed k — /i shadowed variates 
were also derived, while, the moments of this model were presented in [0. In the sequel, the 
model presented in ||T| is used to develop the LCR and AFD equations proposed here. It is worth 
highlighting that the models proposed in Q]| and f2| are related by a simple scaling factor applied 
to the dominant signal component and thus either could be used to arrive at the second-order 
equations presented here. 

The received signal envelope, R, of the n — ft shadowed fading model may be expressed in 
terms of the in-phase and quadrature components of the fading signal such that QJ 
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fl 2 = £(* + &<) 2 + ( y i + ta ) 2 


2=1 


( 1 ) 


where /i is the number of multipath clusters, which is initially assumed to be a natural numbc 




X, and Xj are mutually independent Gaussian random processes with mean E[X,] = E[Yi\ = 0 
and variance E[X 2 ] = E\ Y 2 = a 2 (i.e. the power of the scattered waves in each of the clusters). 
Here pi and q t are the mean values of the in-phase and quadrature phase components of multipath 
cluster i and d 2 = ]X=, p 2 + q 2 . In this model, all of the dominant components are subject to 
the same common shadowing fluctuation, £, which is a Nakagami-m random variable with the 
shaping parameter m used to control the amount of shadowing experienced by the dominant 
components and E[l £ 2 ] = 1. As with the k — p model [[3]|, k > 0 is simply the ratio of the total 
power of the dominant components ( d 2 ) to the total power of the scattered waves ( 2 per 2 ) and the 
mean power is given by E[R 2 ] = r 2 = d 2 + 2pa 2 . While the PDF of R, fn{r), could be obtained 
from |[ 2 l eq. ( 8 )] by expressing the mean power of the dominant component (Q) in terms of n 
and r, that is Q = nr 2 /(n + 1) , for the purposes of this derivation it is obtained from the PDF of 
the instantaneous signal-to-noise ratio ( 7 ) given in JH eq. (4)] via a transformation of variables 
(7 = r 2 7 /r 2 ) as 



( 2 ) 


where FlV) is the gamma function and 1 F } (•: •; •) is the confluent hypergeometric function [fl2ll . 
In this model, m is allowed to take any value in the range m > 0 where m — 0 corresponds 
to complete shadowing of the resultant dominant component and m —> 00 corresponds to no 
shadowing of the resultant dominant component. Of course when m = 00 , the PDF given in (2) 
becomes equivalent to the k — p PDF given in OJ, whereas when rn — 0 and hence k = 0, the 
PDF given in (2) reduces to the Nakagami PDF 0. 

1 Note, this restriction is later relaxed by allowing /r to assume any positive real value. 
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III. Level Crossing Rate 

The level crossing rate of a fading signal envelope, Nr(t) is defined as the expected number 
of times that the envelope crosses a given signal level in a positive (or negative) direction per 
second and is given by [fl3l 


poo 

N R {r)= ff RR (r,f)df (3) 

Jo 

where f is the time derivative of r and f RR {r,r) is the joint probability density of R and /?. If 
we initially hold the shadowing fluctuation constant, the variation of the signal envelope would 
follow a k —n distribution 0. In this instance, from (1), it is easy to see that the k — fi signal 
power can be obtained as the sum of // squared Rice variates. This is a simple but important 
observation as it allows us to show that the PDF of the time derivative of R, denoted as R, 
is zero-mean Gaussian distributed. Letting Z represent a Rice distributed random variable, it 
follows that 


R 2 = J2 Z i- ( 4 ) 

i= 1 

Differentiating both sides of (4) with respect to time we find that 

a = Sk£i£i. (5) 

Knowing that for the Rice case, Z is a zero mean Gaussian distributed random variable with 
variance df = 2" 2 f 2 n a 2 [ 14, eq. (2.104)], where f m is the maximum Doppler frequency, it is 
straightforward to show that in fact & R = & 2 Z . From [3], a 2 = r 2 /2/i(l + k) and therefore 


7 r 2 f 2 r 

n J m' 

At(l + k ) ’ 


Most importantly though, from (5), we can see that R is obtained as a linear transformation of 
Z and thus it can be deduced in a similar fashion to lfT3ll that the PDF of the rate of change of 




6 


the envelope R is uncorrelated with R and thus f r.fSt 1 0 — /rt(' r ) x fh^r)- 

Now considering the shadowed fluctuation of the dominant component separately, which in 
this model is assumed to follow a Nakagami-m distribution. Using the model given in [fl5ll it has 
already been shown that the slope is zero mean Gaussian distributed and its PDF independent 
of the envelope and thus f RR (r,r) = f R (r ) x f R (r). Knowing that the Nakagami-m distribution 
appears as a special case of the k — p distribution, the variance of the slope can also be obtained 
by letting n = 0 in ( 6 ), p — m and interchanging r 2 with D such that a R = n 2 As above, 

to remove the dependency of the formulations on the mean power of the dominant component, 
we substitute Q = nf 2 /{ 1 + k) which gives 


.2 = n 2 flKr 2 
R m( 1 + k ) 


(7) 


Having demonstrated that for both the multipath and shadowing, the variation of the fading 
components are independent of the PDF of the time derivative of the envelope, it now becomes 
possible to rewrite (3) as 


roo 

N R (r) = f R (r ) / ff fi {r)dr ( 8 ) 

J 0 

where f R (r) is the k — ji shadowed PDF given in (2) and f R (r) is the PDF of the rate of 
change of the envelope R. Following the approach taken in [|7]| it seems reasonable to assume 
that the PDF of R is the result of two correlated zero-mean Gaussian random processes. Letting 
R = A + B where A is the rate of change of the envelope due to the multipath component and 
B is the rate of change of the envelope due to the shadowed dominant component, the joint 
density of A and B is given by llT6l 


W*’ i) = SWP7 XHp 


1 /a 2 2 pab b 2 N 

2 ( 1 -P 2 )\^A va°b + v 2 B/ 


\p\ < 1 (9) 


where a a and a B are the variances of the two random variables A and B, and p is the correlation 
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between them. Substituting A = R — B into (9), the integral f r(t) — f-oo f a,b(t — b, b)db can 
be evaluated as m 


fair) = 


[27t(1 - p 2 )(&\ + 2 pa A d B + v 2 b)] 112 


x exp 


f 2 

( O’K 1 - P 2 ) +4p(T j4 Cr B \ 

1 

LO 

l 

to 

^to 

( <7a + 2 pa A d B + cx 2 b j _ 


( 10 ) 


and therefore 



\J( 1 - P 2 ) (b~A + ZP^AO-b + V B )°A 

(dyi(l - p 2 ) + 4 :P&b) 


( 11 ) 


Substituting (2), (6), (7) and (lljjinto (8) and performing the necessary mathematical operations, 
we obtain the LCR of the k — p shadowed fading envelope (normalized to the maximum Doppler 
frequency, f m ) as 


N R (r ) \J 27t (1 — p 2 )p 2 m m ( 1 + /r) M 2 (m + pn + 2 p^/pKmfj 2 

,fJ 


frr 


r (ylf) (ylfft + m) m (i/m (1 — yO 2 ) + 


( 12 ) 


x exp 


p (1 + k) r 2 
^2 


i-Pi rn\ p; 


p?K (1 + k) (r 


pn-\- m \ r / 

For the case when the slopes of the multipath and shadowed components of the received signal 
are uncorrelated (i.e. p — 0), (12) can be further reduced to 


N R (r ) y/2npP 2 m m 2(1 + k) m 2 (m + pn) 2 fr\ 2 ^ 1 


frr 


T (yfi) (/ 1 ft + m)' 


r / 


x exp 


p (1 + k) r 2 


1 F 1 m; /i; 


/r 2 ft (1 + ft) fr \ 


pn-\-m \r J 


- 


(13) 


Fig. 1 shows the shape of the normalized k — p shadowed LCR given in (12) for increasing 
values of m (continuous lines) and decreasing values of p (dashed lines). It is quite evident that 


2 (6) is used in place of while (7) is used in place of a% 
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as the amount of shadowing of the resultant dominant component decreases, i.e. m gets larger, 
the signal crosses lower levels at lower rates. Furthermore, the impact of increasing correlation 
between the slope of the shadowed dominant and multipath signals also acts to cause the signal 
to cross lower levels at lower rates. Fig. 2 shows the normalized LCR of the k — fi shadowed 
fading signal for the special cases when it coincides with the normalized LCRs of the Nakagami 
ffT5l . Rice Ifl4l and k — n lUTTl fading models i.e. K — m — 0 for Nakagami, fi = 1 and m = oo 
for Rice and m = oo for k — ji. 


IV. Average Fade Duration 


The average fade duration (AFD) of a fading signal envelope, T ri (r), is defined as the average 
length of time that the signal spends below the threshold level R and is related to the LCR 
through the relationship |fT5il 


F R (r) 

N R (r)' 


(14) 


As we can see, to calculate the AFD, it is necessary to have an expression for the cumulative 
distribution function, F R (r), of the n—fi shadowed fading signal. As the CDF of the instantaneous 
SNR in k — /i shadowed fading channels has conveniently been derived in [1, eq. (6)], to obtain 
the CDF of the received signal envelope, the same quadratic transformation used to yield the 
PDF of the fading signal is employed which gives 



(15) 


where < T> 2 (*, •; •; •, •) is the bivariate confluent hypergeometric function. Now, the normalized 
AFD of a k — n shadowed fading signal can be straightforwardly obtained by substituting (12) 
and (15) into (14) which gives 






T R (r)f m = 


(v^:(l-p 2 )+4pV^)(l+ K )4 / | i. Ul+Mr 2 n(l+n)m r 2 \ 

. — -— ~ 1 2 I — i, = 2 —, i UK+m \ ^2 ) 

^/27r(l—p 2 )/i2 (m+/ire+2p v //ire?7i) ^ 


exp 


^ p(l+K)r 2 ^ 


_P, ( m - n- M 2 «(l+«) r 2 
4 \ > r"' 1 fin+rri r 2 


(16) 


Again for the case when the slopes of the multipath and shadowed components of the received 
signal are uncorrelated (i.e. p = 0), (16) can be further reduced to 


T R (r)f m 


\/™( 1 +k)^ ( , i. 

—-T-—■—"T=tt>2 ~ m, m,p+ i, 

V27rp^(m+pK) 5 ' ' 


exp 




/x(l+«:)r 2 ^( 1 +Av)mr 2 \ 

f 2 ’ (/i/c+m) f 2 / 

;U 2 /t(l + ^) r 2 \ 
f 2 y 


(17) 


Fig. 3 shows the shape of the k — /i shadowed AFD for increasing values of m (continuous 
lines) and decreasing values of p (dashed lines). As the shadowing of the dominant component 
decreases, the fading envelope spends more time at lower threshold levels. This is in direct 
contrast to the correlation between the time derivative of the multipath and shadowed dominant 
components. In this instance, as the correlation increases , the signal can spend more time at 
lower threshold levels. 


V. A Comparison with Measured Shadowed Fading Channels 

To illustrate the utility of the new equations for modeling shadowed fading channels, they were 
compared with data obtained from two different sets of field measurements which considered 
channels which are known to susceptible to shadowed fading. The first set of measurements 
considered cellular device-to-device communications channels operating at 868 MHz in an 
outdoor urban environment. Full details of the experimental setup are available in [O, [|T8ll . 
This particular scenario considered two persons spaced 10 m apart in an open space between 
three buildings in a built up residential area in the suburbs of Belfast in the United Kingdom. 
Both persons were initially stationary, in direct LOS and had the hypothetical user equipment 
(UE) positioned at their heads. They were then instructed to move around randomly within a 
circle of radius 1 m from their starting points while imitating a voice call. The first person’s 
UE was configured to transmit data packets at a power level of 0 dBm with a period of 70 ms 
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for 80 s while the second person’s UE recorded the received signal strength indicator (RSSI) of 
the received packets. The second set of field measurements considered on-body communications 
channels operating at 2.45 GHz within a highly reverberant environment. Full details of the 
experimental setup are available in llT9ll . The on-body link spanned the left-waist to right-knee 
positions and the measurements were made when the person performed walking on the spot 
movements. In this instance, the complex 5 ' 2 i was sampled with a period of 5 ms for an interval 
of 30 s. 

Fig. 4 shows the empirical level crossing rates for both channels compared to the new equation 
given in (12). All parameter estimates for the k — /i shadowed fading model were obtained using 
the lsqnonlin function available in the Optimization toolbox of MATFAB along with the 
PDF given in (2). It should be noted that both sets of data were normalized to their respective 
rms signal levels prior to parameter estimation. Using these parameter estimates, the maximum 
Doppler frequency and correlation were then obtained by minimizing the sum of the squared error 
between the empirical and theoretical FCR plots. As we can quite clearly see, the normalized 
FCR of the k —[i shadowed fading model provides an excellent fit to the on-body data and a very 
good approximation of the device-to-device channel. To allow the reader to reproduce these plots, 
parameter estimates for both measurement scenarios are given in Table I. In both applications, it 
is quite clear that the fading channel is subject to heavy shadowing with m < 0.55. Furthermore, 
for the device-to-device fading channel, it is apparent that the correlation between the slope of 
the shadowed dominant and multipath components is non-zero. For completeness, Fig. 5 shows 
both the empirical and theoretical AFD of both types of fading channel. Again, the theoretical 
AFD provides an excellent representation of the measured data for the on-body fading channel 
and a good fit for the device-to-device channel for signal levels above -10 dB threshold level. 

VI. Conclusion 

Closed-form expressions for the FCR and AFD of the recently proposed k—h shadowed fading 
model have been presented. These new, very general equations will find use in a wide variety 
of existing and emerging communications applications in which the received signal is subject to 
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shadowed fading such as device-to-device communications, body centric communications and 
land mobile satellite communications. The analytical expressions have been validated through 
reduction to known special cases. It was found that decreasing shadowing of the resultant 
dominant component reduces crossings at low signal levels but at the same time may increase 
fade durations at these levels. Finally the utility of the new formulations has been proven through 
comparison with empirical data obtained for cellular device-to-device and body centric fading 
channels. It has been shown that the second-order statistics of the shadowed fading in both types 
of channel can be adequately described by the equations proposed here. 
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TABLE 1 

Estimated Parameters for Measured Shadowed Fading Channels 


Fading Channel 

k 

P 

r 

m 

fm 

P 

D2D 

1.39 

1.78 

1.14 

0.55 

2.40 

0.29 

On-Body 

0.66 

1.39 

1.03 

0.36 

4.68 

0.05 





r (dB) 


Fig. 1. Normalized level crossing rate for the k — p shadowed fading model with decreasing shadowing of the resultant 
dominant component (continuous lines, p = 0) and with increasing values of the correlation coefficient (dashed lines, m = 1). 
It should be noted that for all of the plots, k = 0.5, ll = 2 and f = 1. 
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Fig. 2. Normalized level crossing rate for the special cases when the k — p shadowed fading model (continuous lines) coincides 
with the normalized LCR of the Nakagami GD (triangle shapes), Rice G3 (square shapes) and k — p G3 (circle shapes) 
models. It should be noted that for all of the plots f = 1 and p — 0 for the k — p, fading model 





-40 -30 -20 -10 0 10 


r(dB) 


Fig. 3. Average fade duration for the k — /j, shadowed fading model with increasing shadowing of the shadowed dominant 
component (continuous lines, p = 0) and with increasing values of the correlation coefficient (dashed lines, m = 1). It should 
be noted that for all of the plots, k = 0.5, p = 2 and f = 1. 
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Fig. 4. Comparison of the theoretical (continuous lines) and empirical (shapes) normalized LCRs for the device-to-device and 
on-body fading channels. All parameter estimates for the theoretical plots are given in Table I. 
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Fig. 5. Comparison of the theoretical (continuous lines) and empirical (shapes) normalized AFDs for the device-to-device and 
on-body fading channels. All parameter estimates for the theoretical plots are given in Table I. 




